Turbulent dynamo action in the high— conductivity hmit: 
a hidden dynamo 

K.-H. Radler and U. Geppert^ 

Abstract. The paper deals with a simple spherical mean-field dynamo 
model of a^-type in the case of high electrical conductivity. A spher- 
ically symmetric distribution of turbulent motions is assumed inside a 
spherical fiuid body surrounded by free space, and some complete form 
of the corresponding mean electromotive force is taken into account. For 
a turbulence lacking reflectional symmetry finite growth rates of magnetic 
fields prove to be possible even in the limit of perfect conductivity, that is, 
the model corresponds to a fast dynamo. In accordance with the theorem 
by Bondi and Gold the dynamo-generated field is completely confined in 
the fluid body. 



1. Introduction 

According to a finding by Bondi and Gold (1950) the magnetic multipole mo- 
ments which result from electric currents in a perfectly conducting fluid occu- 
pying a simply connected body cannot grow boundlessly. They remain for any 
fluid motion within bounds determined by the initial magnetic flux through the 
surface of the body. In a simple spherical mean-field dynamo model of a^-type 
proposed by Krause and Steenbeck (1967), however, the magnetic field grows 
endlessly both inside and outside the fiuid body even in the limit of perfect 
conductivity. This confiict has been resolved in a more sophisticated model by 
Radler (1982) that properly considers the specific structure of the mean electro- 
motive force resulting from the constraints on the fluid motion at the boundary 
of the body. It was shown analytically that this structure indeed ensures the 
boundedness of the magnetic multipole moments and thus excludes an infinite 
growth of the magnetic field outside the fluid body in the high-conductivity 
limit. However, the question remained open whether or not a dynamo may 
work inside the fluid body, which then would have to be invisible, or hidden, in 
the sense that its magnetic held is completely confined inside this body. Only 
a few arguments were given which support the conjecture that such a dynamo 
may exist. 

In between the issue of the adjustment of a dynamo in the high-conductivity 
limit to the constraints posed by the Bondi and Gold theorem was also ad- 
dressed in papers by Hollerbach, Galloway and Proctor (1995 and 1998). They 
considered a dynamo in a spherical fiuid shell with a chaotic flow surrounded 
by insulating space inside and outside and demonstrated that it works also in 
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the limit of high conductivity, in which then the magnetic field in the insulating 
spaces vanishes. 

In the present paper we return to the mean-field model proposed by Radler 
(1982). After a few short explanations concerning the case of high conductivity 
and the mean-field concept in dynamo theory (Sections 2 and 3) and the sim- 
ple model by Krause and Steenbeck (1967) (Section 4) we deliver a systematic 
treatment of our model and present numerical results giving evidence for the 
existence of a hidden fast dynamo in the above sense (Section 5). Finally we 
discuss a few conclusions (Section 6). 

2. Dynamo action in the high conductivity limit 

We consider here dynamo models consisting of an electrically conducting fluid 
occupying a simply connected region surrounded by free space. Let us assume 
that the magnetic flux density, B, is governed by the induction equation 

rjV^B + 'V X {ux B) -dB/dt = 0, V-B = 0, (1) 

inside the fluid body, continues as a potential field 

B = -V(p, A(j) = 0, (2) 

in outer space and vanishes at infinity. As usual, u means the velocity of the 
fluid and r] its magnetic diffusivity. 

For steady u we may expect solutions B = ^(B exp(pt)), where S is a 
steady complex field and p a complex quantity independent of space and time 
coordinates. The real part pr of p is the growth rate of B. We speak of a 
dynamo if there is at least one solution B with a non-negative pr- 

For many applications to cosmic objects the limit of high electrical con- 
ductivity, that is 7/ — > 0, deserves special interest. To define this limit more 
precisely, let us measure all length in units of L being a typical scale of it or B, 
and the time in units of L/U, with U being a typical magnitude of u. Then (|^) 
applies with replaced by -Rm"^, where i?m is the magnetic Reynolds number 
defined by 

Rm = UL/7], (3) 

and the high-conductivity limit corresponds to Rm oo. 

We mention two important aspects of dynamo action in the high-conduc- 
tivity limit. To explain the first one we compare dynamos which differ only in 
the value of Rui- If then the largest growth rate pr remains positive and takes 
a finite positive value as i?m — > oo we speak of a fast dynamo, otherwise, that 
is, if Pr tends to zero or to a negative value, of a slow dynamo. 

The second aspect is that in the case of perfect conductivity, Rm^^ = 0, the 
magnetic field has to satisfy the requirements posed by the theorem by Bondi 
and Gold (1950). For the sake of simplicity we suppose the fluid body to be 
a sphere and formulate this theorem as proposed by Radler (1982). We use 
spherical coordinates r, -i?, ip and represent the potential (p introduced with (^) 
in the form 

oo n 

= E E Cr-("+^)yr(^,'/^) (4) 

n=l m=—n 
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with spherical harmonics defined by Y^{'d, ip) = P^{cos'd) expiiraip) where 
the are associated Legendre polynomials. The are complex constants 
satisfying c~"^ = c^* which define multipole moments, e.g., those with n = 1 
the dipole moment, with n = 2 the quadrupole moment etc. Depending on the 
motions at the boundary of the fluid body the may vary in time. According 
to the theorem by Bondi and Gold, however, they are bounded in the sense that 

ICI<C, (5) 

with given by the initial distribution of the magnetic flux at the boundary. 
This in particular excludes any exponential growth of the . 

A third aspect which concerns the conditions to be satisfied by B at the 
boundary will be discussed below. 



3. The mean— field approach 

In cases in which the magnetic field and the fluid motion are of turbulent nature, 
or show by other reasons complex structures in space or time, the mean-field 
approach to dynamo models has proved to be useful; see, e.g., Krause and 
Radler (1980). In this approach both the magnetic flux density B and the 
fluid velocity u are considered as a sum of a mean field, B or tl, defined by a 
proper averaging procedure, and a fiuctuating part, B' or u' . Provided that the 
Reynolds averaging rules apply, we may conclude from equations (|l|) that 

riV'^B + V ■K{uxB + e)-dB/dt = Q, V-B = 0, (6) 

inside the fiuid body, where £ is an electromotive force due to fluctuations, 

£ = u' xB' . (7) 

For a given motion ^ is a linear functional of B. Under the usually accepted 
assumption of sufficiently weak variation of B in space and time it can be rep- 
resented as 

£i = ttijBj - bijkdBj/dxk . (8) 

Here we rely on Cartesian coordinates and use the summation convention. The 
tensors aij and bijk are, apart from rj, determined by u and u'. 

If, for instance, u = and u' represents a homogeneous isotropic turbulence 
we have 

e = aB-l3V xB (9) 

with coefficients a and f5 determined by u' or, what is the same here, by u and 
being independent of space coordinates. The first term on the right-hand side 
describes the o-effect, that is, the occurrence of a mean electromotive force par- 
allel or antiparallel to B but vanishes if the turbulence is reflectionally symmet- 
ric. The second term gives rise to introduce a mean-field conductivity different 
from the usual one. Several results are available concerning the connection of 
a and f3 with the properties of the it-field. We mention in particular those 
obtained in the second-order correlation approximation; see Krause and Radler 
(1980). In the high-conductivity limit, in this context defined by rjTc/X'^ — > 
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with Ac and Tc being correlation length and time of the u-field, a and (3 then 
take the values a^^^ and P^^^ given by 

a^^> = \ u{x,t) ■ {V xu{x,t-T))dT 
Jo 

(10) 

roo 

/3(0) = i / u{x,t)-u{x,t-T)dT. 
Jo 

In general both a and /3 do not vanish in this limit. 



4. A (too) simple model 

One of the simplest mean-field dynamo models, which can be treated analyt- 
ically, has been proposed by Krause and Steenbeck (1967). The fluid body is 
supposed to be a sphere of radius R, again surrounded by free space. Any mean 
motion of the fluid is ignored, u = 0, and thinking of homogeneous isotropic 
turbulence the electromotive force £ is taken in the form (|9|) with constant a 
and p. Specifying now equations (|6|) to this case we write simply B instead 
of B and use R and /{rj + j3) as units of length and time. Looking then for 
solutions B varying like exp(pt) with t we may reduce these equations to 

V^B + CV X B -pB = Q, V-S = 0, (11) 

where C is a dimensionless measure of the a-effect, 

C = aR/{'n + [3). (12) 

In this model equation (^) was completed by the requirements that B is con- 
tinuous across the boundary and has the structure given by @ and ^ in outer 
space. The problem posed in this way has been investigated for axisymmetric 
S-fields by Krause and Steenbeck (1967) for the steady case and by Voigtmann 
(1968) for the time-dependent case. A detailed treatment of the general case in- 
cluding non-axisymmetric time-dependent fields is given in Krause and Radler 
(1980). We only mention here a few particular results. There are independent 
solutions B which possess the form of single multipole fields, that is dipole fields, 
quadrupole fields etc. in outer space. In all cases p is real, it takes a negative 
value for C = and increases monotonously with | C | , runs through zero for 
some marginal value of | C | and behaves like as C oo. Comparing all 
these solutions we find the smallest marginal value for a solution of dipole type. 
Denoting this value by Ccrit we have 

Ccrit = 4.493 . (13) 

Consider now the high-conductivity limit. Then the condition of dynamo 
action resulting from (12) and (13) takes the form \ a.^^^ \R/f5^^'^ > 4.493. It 



seems well possible to satisfy this condition, that is, the dynamo may well work 
in this limit. The dimensional growth rate is then given bmce 
in the case of a dynamo with growing magnetic fields p must be positive, this 
growth rate is positive too, that is, we have a fast dynamo. 
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The structure of the magnetic field depends, if its multipole character is 
specified, only on the value of the parameter C irrespective of that of ??. The 
field is non-zero both inside and outside the fluid body, and this cannot change 
as 7? ^ 0. If a dynamo works and the field grows exponentially in time then 
the c™ introduced with do so, too. This, however, is in conflict with the 
Bondi-Gold theorem. 

5. The more sophisticated model 
5.1. Basic equations 

One of the shortcomings of the model considered so far, which might be the rea- 
son for this conflict, is the assumption of an electromotive force £ in the form (^) 
with constant a and /3, which can be justified for homogeneous isotropic turbu- 
lence only. Near the boundary of the fluid body the turbulence must necessarily 
deviate from homogeneity and isotropy, and £ has to take a more complex form. 
As already mentioned, in a paper by Radler (1982) a modified model was studied 
with a spherically symmetric distribution of the turbulence. Then the radial di- 
rection necessarily occurs as preferred direction in the turbulence and therefore 
we have a more complex form of £. It was shown that with this modification 
the conflict is resolved. More precisely, with these assumptions any growth of 
the magnetic field at the boundary of the fluid body or in outer space can be 
ruled out. 

Turning now to this model we assume again that the mean magnetic flux 
density B is governed by equations (^) in a spherical fluid body r < R and 
continues as a potential field as given by (P), or (Q), in the outer space r > R. 
We further exclude any mean motion, u = 0, and specify the electromotive 
force S in accordance with a spherically symmetric turbulence. For the sake of 
simplicity we write in the following again B instead of B, and u instead of u', 
and we rely again on spherical coordinates r, -d, ip. 

Spherical symmetry of the turbulence is understood here in the sense that 
all averaged quantities determined by the it-field are invariant under arbitrary 
rotations of this field about arbitrary axes through the center r = of the fluid 
body. From this definition we may conclude by standard reasoning (see, e.g., 
Krause and Radler (1980)) that the electromotive force S must have the form 

E = — a\B — a2{r ■ B)r — -yf x B 

-piV X B- (32{f- (V X B))f-5f X (V X S) (14) 
- Pl{f ■ (VBY) - /35(f • (f • {VBy))f - S'f X (f • {VBY) 

with scalar coefficients oi, 02, (5^ determined by u and depending on position 
only through r but not through {} ot ip. Here f is the radial unit vector and 
(VBy the symmetric part of the gradient tensor of B so that, if we refer to 
Cartesian coordinates, (f • (VBY)- = \fj{dBi/dxj + dBj/dxi). Some more 
details concerning the derivation of (jl4) are given in Radler (1982). If the 
turbulence is reflectionally symmetric about planes containing the center of the 
body ai,a2,S, (31 and are equal to zero. In the case of homogeneous isotropic 
turbulence only the coefficients ai and Pi can be non-zero and the others have 
to vanish, so that we return to (P). 
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Later we will also assume that the turbulence is steady, that is, that all 
averaged quantities depending on u are invariant under shifts along the time 
axis. In this case the coefficients ai,a2, ■■■,S^ are independent of time. 

Even if we accept again the second-order correlation approximation the 
calculation of the coefficients ai,a2, ...,5'^' for points near the boundary of the 
fluid body is, at least for finite conductivity, rather complex. We note here only 
results for the high-conductivity limit, defined as above by rjTc/X'^ 0, which 
were already obtained in a slightly different form by Radler (1982). Denoting 

the mentioned coefficients in this limit by af\4°\...,5'-W we have 



= Ub± + bii), /3f = -i<5'^(°) = Ub± -h), (15) 



and 



ay = J 'U||(a;,t) • (V X tt(£C, i — r))|| dr , 

/•oc 

ail = / u\\{x,t — t) ■ CV X u(x.t))u dr , 
Jo " 

POO 

aj_ = i / u±{x, t) • (V X u{x, t — r))_|_ dr , 
Jo 



bn = / uu{x.,t) ■Uu(x,t-T)dT, (16) 

Jo 

/■oo 

b± = h u±{x,t) ■ u±{x,t — t) dr , 
Jo 

c = j (^{f ■ u{x,t)){'V ■ u{x,t - t)) - {f ■ u{x,t - t)){V ■ u{x,t))^ dr , 

/•oo 

d= f ■ {u{x,t) X u{x,t — t)) dr , 
Jo 

with = (f • u)f and u± = it — and analogous definitions of (V x t/)| 

and (V X u)j_. According to our assumptions the a^^\a2'\ ■■■■,S^^^\ like the 
averaged quantities under the integrals, depend on x via r only. 

It seems natural to assume that fty and b± are non-negative everywhere. At 
the boundary of the fluid body we have f • tt = and therefore 

ay = ay = 6|| = db\\/dr = c = at r = 1 . (17) 

5.2. Reduction of the basic equations 

Let us now represent the magnetic flux density B as a sum of a poloidal and 
toroidal part, 

S = - V X (r X VS) - {r/R) x VT , (18) 



where r = rf, and expand the defining scalars S and T in series of spherical 
harmonics ip). It can easily be followed up that the equations and condi- 

tions governing B imply no coupling between contributions to B differing in n 
or m so that we may restrict ourselves to the simple solutions defined by 

S = S{r, V),T = T{r, , n > 1 , \m\ < n . (19) 



Due to the factor R in ( |18D the dimensions of S and T coincide. Preparing 
the definition of dimensionless quantities we introduce the constants q'^ and /J" 
with the dimension of a velocity and a magnetic diffusivity. We will assume 
homogeneous isotropic turbulence in some central region of the fluid body and 
identify these constants with a^^^ and P^^^ given by (0) for this region. In 
the following we measure all lengths in units of R, the time in units of R? / (3^ , 
further the ai, 02, ay, dy, a_L in units of , the 5,/?[,/32,'i in units of a^R, the 
/^i) /32 5 (^^j ^|[, ^± in units of and 7,c in units of (3^ /R. 

Using standard methods we may then reduce equations (^) for B to 



eDS + Us + CUt - dS/dt = 
eDT + CVs + Vt- dT/dt = 



for r < 1 (20) 



with 

£ = T^Z/jO , C = a°R/p° , (21) 



^ dHrf) n(n+l) ^ ^^2) 



and 



Us = l'-^ + PiDS - 1,(23 - r^^ - nin + 1) S) 



Ut = -air+-^^ + ^(2T-^^' 
r dr 2r dr 



(23) 



Idf d{rS)\ .n(n + l) 

Vs = — - (ai + 02) 2 ^ ~ -^i^rDS) 

r or \ or / r or 



Idi^ Id d{rT) 

= ^ + — ) 

r or r or or 

n{n+l) 1 d d{rT) 
-(/3i+/32)— :^T---(5 (2r-^— )). 



(24) 
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Equations (|2|) for the outer space are equivalent with 



5 = 5(r = l)/r"+S T = for r>l. (25) 

Note that S{r = 1) coincides with c^/K^~^^ for the chosen n and m. 

As a consequence of the divergence relation for B its normal component and 
thus S has to be continuous across the boundary r = 1. For finite conductivity 
of the fluid, that is e > 0, we may exclude surface currents so that the tangential 
components have to be continuous too, and so dS/dr and T. Together with ( p5| ) 
we may conclude that 

dS/dr + (n + 1)5 = T = at r = 1 . (26) 

The original problem of the determination of B occurs then as the problem 
of solving the equations ([20|), completed by ( p3|) and (p^), with the boundary 



conditions (26) 



In the case of perfect conductivity, e = 0, again the continuity of S has 
to be required. However, surface currents can no longer be excluded, which 
correspond to discontinuities of dS/dr and T. More precisely, such currents are 
given by {1/ iiR^){[dS/dr]f x Vi^ + [T]Vy^), where ji means the magnetic 
permeability of the fluid and [/] stands for /(r = 1 + 0) — /(r = 1 — 0). In 
this case the equations (^0|) have to be completed by conditions which we will 
formulate later. 

5.3. Compatibility with the Bondi Gold theorem 

Remaining with the case of infinite conductivity we first demonstrate, repeating 
the ideas described by Radler (1982), that our model is compatible with the 
Bondi-Gold theorem. For this purpose we consider the first equation (20), com- 



pleted by (|23D, for r — s- 1. Putting there e = 0, inserting ai,a2, ■■■■,8'' according 
to (|l5|) and using ([l7|) we find 



dS/dt = -n{n + 1) b^S at r = 1 , (27) 

and hence 

S{r = 1, t) = S{r = 1, 0) exp ( - \n{n + l)h^{r = l)t) . (28) 

Remembering that S{r = 1) = c™/i?"^^ and that is non-negative we see 
that (|5|) is indeed satisfied. 

5.4. Further specification of the model 

We now turn our attention to the question whether in the limit of infinite con- 
ductivity magnetic fields can grow inside the fluid body. We rely on equations 
(|0|) to dH), put e = and insert ai, 02, • • • , 5'' according to (^). With the idea 



to have a model with a minimum of free parameters we choose simply 

an = 5|| = 6|| = / , a_L = &± = 1 , c = d = (29) 



with 



. , 1 , for < r < 1 - g 



This corresponds to the assumption of homogeneous isotropic turbulence inside 
the spherical region < r < 1 — q and of inhomogeneous anisotropic turbulence 
in the surrounding shell 1 — q < r < 1. Note that / = df /dr = d'^f/dr'^ = at 
r = 1 and therefore (17) is satisfied. Finally we assume that there is initially no 
magnetic field at the boundary of the fluid body, that isS' = r = Oatr = l. 
Due to ( p7[) w e have then S = at r = 1 at any time. Starting from the second 
equation (pO|) together with (p^ ) and using ([T5|), p^ ) and (30) we arrive at a 
relation for T analogous to (|27| ) from which we may conclude that also T = 
at r = 1 at any time. The behaviour of B is then governed by (|20|) , completed 
by (p^, (p^), (p^), ( p9D and ([30|), and the role of boundary conditions is taken 
byS' = r = Oatr = l. (We note that the remarks concerning the boundary 
condition for S in Radler (1982), Section 3.6, are incorrect.) 



5.5. Numerical results 

The problem posed in this way has been investigated numerically. There are 
indeed solutions B, or S and T, which for sufficiently large | C \ grow exponen- 
tially in time. The growth rate p introduced as above again proves to be real. 
The fastest-growing solution for a given | C | is always one of dipole type, that 
is n = 1. We restrict our attention in this paper to such solutions and, further- 
more, to the case q = 0.2. The dependence of the growth rate p on [ C | is shown 
in Figure |^. Clearly dynamo action is possible if | C | > Ccrit with 



Ccrit = 5.002 . 



(31) 
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Figure 1. The growth rate p versus | C | 

Let us define an effective radius i^eff of our model by equating of a°i?efj//3° 
to the value of Ccrit for the simple model discussed above given by (|13|), that 
is a^Res/f3^ = 4.493. Comparing this with the corresponding relation for the 
modified model under consideration, that is a R/pO = 5.002, we find Res/R = 
0.90. This means that the shell with the fractional radius between 0.8 and 1 
showing deviations from a homogeneous isotropic turbulence is less efi^ective for 
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Figure 2. The energy densities and of the poloidal and toroidal 
parts of the magnetic field in arbitrary units in dependence of the 
fractional radius r. Left panel C = 8, right panel C = 12. 

dynamo action than the central region of the fluid body without such deviations. 
Figures |2| and |3| show the radial distribution of the energy densities of the poloidal 
and the toroidal parts of the magnetic field and field pictures for two examples 

with C > Ccrit- 




Figure 3. Field lines of the poloidal part (right) and isolines of the 
toroidal part (left) of the magnetic field. The solid isolines correspond 
to positive, the broken ones to negative values of the (/j-component of 
the field. Left panel C = 8, right panel C = 12. 



10 



It would be interesting to study a sequence of dynamo models with positive 
£ approaching zero. In this context we must have in mind that e occurs not only 
with the dissipation terms eDS and eDT in (|20| ). In addition the coefficients 
«!, 02, • ■ ■ 5 entering Us,Ut, Vs and Vr via (23) and depend on the mag- 
netic diffusivity, that is, on a parameter proportional to e. There are, however 
no results available which describe this dependence in the neighbourhood of a 
jump of the magnetic diffusivity as it occurs at the boundary. 

Several numerical calculations have been carried out on the basis of equa- 
tions ( pO| ) including the dissipation terms but using expressions for cti , 0:2 , ■ ■ ■ , i^i 
as given by (p^) and ([T^), that is, for infinite conductivity. Of course, in these 
calculations the boundary conditions (26) have been used. Indeed the solutions 
obtained in this way approach the corresponding ones with e = in a very sat- 
isfying manner as e — > 0. The only difference occurs in the distribution of the 
electric currents near the boundary, what is understandable since for e = a 
part of them occurs as surface currents. 

It is interesting to observe how the magnetic energy concentrates itself more 
and more inside the conducting body as e ^ 0. Table |l] shows for one example 
the dependence of the ratio of the energies -Eout and Etot in the outer and in all 
space on e. 



e Eout/Etot 



10-2 8.6 • 10^^ 
10-3 4.9 • lO-'^ 
10-^ 4.8 • 10-9 



Table 1. The energy ratio i?out/-E'tot in dependence on e for C = 12. 



6. Conclusions 

As already shown in an earlier paper by Radler (1982), the mean-field model 
reconsidered here, in contrast to a more simplified one, meets the requirements 
posed by the Bondi-Gold theorem. It was, however, only conjectured but not 
really demonstrated that it admits dynamo action inside the fluid body in the 
high-conductivity limit. The present paper clearly demonstrates the possibility 
of a fast dynamo inside the fluid body, whose magnetic field is then, as required 
by the Bondi-Gold theorem, completely confined in the fluid body and in that 
sense invisible from outside. 

The model is also in another respect, which is not necessarily connected 
with the high-conductivity limit, of interest for the fundamentals of mean-field 
dynamo theory. In most of the mean-field dynamo models elaborated in view 
of cosmic objects only a few contributions to the mean electromotive force have 
been taken into account, and others were cancelled in the vague hope that they 
are of minor importance. We demonstrate here the possibility of dynamo ac- 
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tion in an idealised model which is consistent in the sense that it includes all 
contributions to the electromotive force which occur under this idealisation. 

It remains to be discussed whether the concentration of the dynamo-ge- 
nerated magnetic field in the interior of a highly conducting body, which was 
demonstrated above, may indeed occur in astrophysical bodies. 
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